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Abstract 

These notes briefly summarize the lectures for the Summer School " Optimal transportation: 
Theory and applications" held by the second author in Grenoble during the week of June 22- 
26, 2009. Their goal is to describe some recent results on Brenier's variational models for 
incompressible Euler equations [1, 2, 4]. 

1 Euler incompressible equations and Arnold geodesies 

Let D denote either a bounded domain of M. d or the d-dimentional torus T d := W 1 /If 1 . We consider 
an incompressible fluid moving inside D with velocity u. The Euler equations for u describe the 
evolution in time of the velocity field, and are given by 

d t u + (u ■ V)« = -Vp in [0, T] x D, 
divw = in[0,T]xD, 

coupled with the boundary condition 

u ■ v = on [0, T] x dD 

when D ^ T d . Here p is the pressure field, and arises as a Lagrange multiplier for the divergence-free 
constraint on the velocity u. 

If u is smooth we can write the above equations in Lagrangian coordinates: let g denote the flow 
map of u, that is 

g(t,a) = u(t,g(t,a)), 
g(0,a) = a. 

By the incompressibility condition, and the classical differential identity 

^-detV a g(t,a) = div u(t,g(t, a)) deW a g(t, a), 
at 

(here and in the sequel div denotes the spatial divergence of a possibly time-dependent vector field) 
we get detV a g(t, a) = 1. This means that g(t,-) : D — > D is a measure-preserving diffcomorphism 
of D: 

g(t,-)#n D =ii D (i.e. f i D (g{t 1 -)- 1 {E))= f i D (E) VE). 
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Here and in the sequel f#[i is the push-forward of a Borel measure fi through a map / : X — > Y (i.e. 
f Y 4> df#[i = J x 4>o f d/i for all Borel bounded functions <fi : Y — > R), and /id is the volume measure 
of D, renormalizcd by a constant so that /j,d(D) = 1. 

Writing Euler's equations in terms of g we obtain an ODE for 1 1-> g(t) in the space SDiff (£>) of 
measure- preserving smooth diffeomorphisms of D: 

g(t, a) = -Vp (t, g(t, a)) (t, a) G [0, T] x D, 

g(0,a)=a a e D, (1) 

g(t, ■) S SDifT(Z)) ie[0,T]. 

1.1 Weak solutions to Euler's equations 

In the case d = 2, existence of distributional solutions can be proved through the vorticity formu- 
lation: setting uit(-) = curlu(£, •), so that u(t, ■) = V- L A _1 w t , the Euler equations can be read as 
follows: 

—uit(x) + div (u>t(x)u(t, x)) = 0. 

Formally, this equation preserves all V norms of solutions, and indeed existence is not hard to 
obtain if ujq 6 L p for 1 < p < oo. Delort improved the existence theory up to L 1 or measure 
initial conditions wo whose positive (or negative) part is absolutely continuous, and it is still open 
the problem of getting a solution for all measure initial data. As shown by Yudovitch [15, 16] 
uniqueness holds for p = oo, while it is still open in all the other cases. 

In the case d > 2 much less is known: no general global existence results of distributional 
solutions is presently available. 

1.2 Arnold's geodesic interpretation 

At least formally, one can view the space SDiff (D) of measure-preserving diffeomorphisms of D as an 
infinite-dimensional manifold with the metric inherited from the embedding in L 2 (D;M. d ), and with 
tangent space made by the divergence-free vector fields. Using this viewpoint, Arnold interpreted 
the ODE (1), and therefore Euler's equations, as a geodesic equation on SDiff (D). Therefore one 
can look for solutions of Euler's equations on [0, 1] x D by minimizing the Action functional 

*(s) ~ I I \\g{t 1 x)\ 2 d^ D (x)dt 

JO JD 1 

among all paths git,-) : [0,1] — > SDiff(D) with (?(0, •) = / and g(l,-) = h prescribed (typically, 
by right invariance, / is taken as the identity map i). Ebin and Marsden proved in [10] that this 
problem has indeed a unique solution when /io/ _1 is sufficiently close, in a strong Sobolev norm, to 
i. We shall denote by 6(f, h) the Arnold distance in SDiff (D) induced by this minimization problem. 

Of course, this variational problem differs from Euler's problem, because the initial and final 
diffeomorphisms, and not the initial velocity, are prescribed. Nevertheless, the investigation of 
this problem leads to difficult and still not completely understood questions (typical of Calculus of 
Variations) namely: 

(a) Necessary and sufficient optimality conditions; 

(b) Regularity of the pressure field; 

(c) Regularity of (relaxed) curves with minimal length. 

Before describing some of the main contributions in this field, let us recall some "negative" results 
that motivate somehow the necessity of relaxed formulations of this minimization problem. 
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1.3 Non-attainment and non-existence results 

Shnirelman [12, 13] found the example of a map g 6 SDifF([0, l] 2 ) which cannot be connected to i 
by a path with finite action, i.e. S(i,g) = +oo. Furthermore, he proved that for h € SDiffQO, l] 3 ) of 
the form 

h(xi, x 2 , x 3 ) = (g 1 (xi,x 2 ), g 2 (xi,x 2 ), x 3 ), with (g x , g 2 ) = 9 as above, 

S(i,h) is not attained, i.e. no minimizing path between i and h exists (although there exist paths 
with a finite action). This fact can be easily explained as follows (see also [ti, Paragraph 1.3]): since 
there is no two dimensional path with finite action connecting i to g while in 3 dimension it is known 
that the minimal action is finite [12], if a minimizing path t i— >• g(t) exists then it has a non-trivial 
third component, i.e. g 3 (t,x) ^ x 3 . Set 77(23) := min{2a;3,2 — 2x3}, and let u denote the velocity 
field associated to g, i.e. u = g o g^ 1 . Then it is easily seen that the velocity field 

!ui {xi, x 2 , ry(x 3 )) 
U 2 {xi, X 2 , T)(x 3 )) 
^U 3 (x 1 ,x 2 ,T](x 3 )) 

induces a path g which still joins i to h, but with strictly less action (since u 3 is not identically 
zero). This contradicts the minimality of g, and proves that there is no minimizing path between i 
and h. (See also [8, Paragraph 1.3].) 

Let us point out that the above argument shows that minimizing sequences exhibit oscillations 
on small scales, and strongly suggest the analysis of weak solutions. 



1.4 Time discretization, minimal projection and optimal transport 

Before describing the concept of relaxed solutions to the Euler equations introduced by Brenier, let us 
first see what happens when one tries to attack the above variational problem by time-discretization: 
assume D C M d , and fix go^i G SDiff(D). We want to find the "midpoint" gi/ 2 between go and gi, 
that is we consider 

Up to rearranging the terms and removing all the quantities independent on g , the above problem 
is equivalent to minimize 

,9o + .91 
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L 2 (D;R<i)' 



i.e. we have to find the L^projcction on SDiff(£>) of the function £ L 2 (D;R d ). Since the set 

SDiff(D) is neither closed nor convex, no classical theory is available to ensure the existence of such 
projection. 

In order to make the problem more treatable, let us close SDiff(Z?): as shown for instance in [9], 
if D = [0, l] d or D = T d then the L 2 -closure of SDiff {D) in L 2 (D; R d ) coincides with the space S(D) 
of measure preserving maps: 

S(D) :={g:D^D: n D [g-\A)) = /x D (A) VA e B(D)} . 

Then the general problem we want to study becomes the following: given h £ L 2 (D;M. d ), solve 

\h-s\ 2 d^i D . (2) 



mm 

ses(D) J D 



As in the classical optimal transport problem, one can consider the following Kantorovich relaxation: 
denoting by n(R d ) the set of probability measures on M. d x K d with first marginal fin and second 
marginal v := h^^io, we minimize 



I re — 2/ 1 d'j(x,y). (3) 
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Assume the non- degeneracy condition v <C dx. Then we can apply the classical theory of optimal 
transport with quadratic cost for the problem of sending v onto \in [6] : there exists a unique optimal 
transport map V<f> : M. d — > M. d such that (V^>)#f = (Id- Moreover the unique optimal measure 7 
which solves (3) is given by 

7 = (V0 x Id)#i/. 

Then it is easily seen that the map 

s := o ft, 

belongs to S'(-D) and uniquely solves (2) (see [G] or [14, Chapter 3] for more details). 

2 Relaxed solutions 

In the last paragraph we have seen how the attempt of attacking Arnold's geodesies problem by 
time discretization leads to study the existence of the L 2 -projection onto SDiff(D), and that the 
projection of a function h onto its closure S(D) exists and is unique whenever h satisfies a non- 
degeneracy condition. Instead of going on with this strategy, we now want to change point of view, 
attacking the problem by a relaxation in "space" . 

Two levels of relaxation can be imagined: the first one is to relax the smoothness and injectivity 
constraints, and this leads to the definition of the space S(D) of measure-preserving maps. However, 
we will see that a second level is necessary, giving up the idea that g(t, •) is a map, but allowing it 
to be a measure preserving plan (roughly speaking, a multivalued map). This leads to the space 

T(D) := {ry £ &{D x D) : rj{A x D) = fi D {A) = rj(D x A) VA e B{D)} . 

The space S(D) "embeds" into T(D) considering 

S(D)3g^(ixg) # fi D eT(D). 

Conversely, any r\ £ F(D) concentrated on a graph is induced by a map g € S(D). 

Even from the Lagrangian viewpoint, it is natural to follow the path of each particle, and to 
relax the smoothness and injectivity constraints, allowing fluid paths to split, forward or backward 
in time. These remarks led in 1989 Brcnier to the following model [5]: let 

fi(£>) := C ([0, l];D) , e t (cu) := t G [0, 1]. 

Then, denoting by &(Q,(D)) the family of probability measures in H(D), we minimize the action 
functional 

<*(n)-= I \l \u\ 2 dtdr,(cj), ijG^(fi(D)) 
Ja(D) L Jo 

with the endpoint and incompressibility constraints 

(e , ex) # J7 = (ix h)#n D , (e t )#TJ = Vt e [0, T\. 

In Brenier's model, a flow is modelled by a random path with some constraints on the expectations of 
this path. As we will see below, this problem can be recast in the optimal transportation framework, 
dealing properly with the incompressibility constraint. 

Classical flows g(t,a) induce generalized ones, with the same kinetic action, via the relation 
V = {$> g )#HD, with 

$ g :D^n(D) 7 <f> g (a):=g(;a). 

In this relaxed model, some obstructions of the original one disappear: for instance, in the 
case D = [0,l] d or D — T d it is always possible to connect any couple of measure preserving 
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diffcomorphism by a path with action less than \fd. Actually, this allows to prove that finite-action 
paths exist in many situation: as shown in [1, Theorem 3.3], given a domain D for which there 
exists a bi-Lipschitz measure-preserving diffcomorphism $ : D — > [0, by considering composition 
of generalized flows with $ one can easily constructs a generalized flow with finite action between 
any /io,/ii € SDiff(D). Moreover, standard compactness/lower semicontinuity arguments in the 
space &(n(D)) provide existence of generalized flows with minimal action. 

2.1 Eulerian-Lagrangian model 

Coming back to the relaxed model described above, we observe that the endpoint constraint (eo, ei)#Tf = 
(i x h)#fi£> cannot be modified to deal with the more general problem of connecting / 6 S(D) to 
h G S(D): indeed, by right invariance, this is clear only if / is invertible (in this case, one looks 
for the optimal connection between i and h o J -1 ). These remarks led to a more general model, 
which allows to connect r\ = rj a <8 fJ-D to 7 = j a ® fj,r> [2]. (Here we are disintegrating both the 
initial and final plan with respect to the first variable.) The idea, which appears first in Brenier's 
Eulerian-Lagrangian model [8] is to "double" the state space, adding to the Eulerian state space 
D a Lagrangian state space A. Even though A could be thought as an identical copy of D, it is 
convenient to denote it by a different symbol. 
Let 

Q*(D) := fl(D) x A. 

Then, consider probability measures r/ = r\ a ® \iu in f2*(D): this means that r\ has ^ as second 
marginal, and that 

<j>(u,a)dri(w,a) = / ( / <f>(cj, a) dr] a ) d^ D (a) 
J J a \Jn(D) / 

for all bounded Borel functions <j> on £l*(D). 
Again, one minimizes the action 

■= I If \u\ 2 dtd V {LJ,a) 
Jn*{D) 1 Jo 

with the incomprcssibility constraint (e t )#r) = \id for all t (here e t (w,a) = u)(t)) and the family of 
endpoint constraints: 

{e )#V a = 7a, {ei)#V a = Va for /i D -a.c a e D. 

As in the previous section, we are using rj a <£> [Id and j a ® M£> to denote the disintegrations of 77 and 
7 respectively. 

Denoting by 8(n,~f) 2 the minimal action, it turns out that one can define natural operations of 
reparameterization, restriction and concatenation in this class of flows. These imply that (6, T(D)) 
is a metric space. 

Indeed, it is proved in [1] that it is complete and a length space, whose convergence is stronger 
than weak convergence in &(D x D). 

2.2 Motivation for the extension to T(D) 

Even for deterministic initial and final data, there exist examples of minimizing geodesies rj that are 
not deterministic in between: this means that (eo, et)#T] € T(Z?) \ S(D), t € (0, 1). 

To show this phenomenon, consider the problem of conncctingup to additive constants in D = 
-Bi(O) C M 2 the identity map i to —i. For convenience, up to a reparameterization, we can choose 
the time interval as [0,7r]. Two classical solutions are 

[0, 7r] 3 t H> (xi cos±t + X2 sin±£, x\ sin±£ + X2 cos±i), 
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corresponding to a clockwise and an anti-clockwise rotation. 

On the other hand, one can consider the family of maps u x ,e connecting x to —x 

u x .e(t) := a; cost + i/l - \x\ 2 (cos9, sin#) sint 9 e (0,7r) (4) 

and define 77 := (u x , e )^^\D x if 1 [(0, 2tt)) . 

It turns out that 77 is optimal as well, and non-deterministic in between. Moreover, as shown 
in [4], it is possible to construct infinitely many other solutions to the above minimization problem 
which are not induced by maps. For instance, one can split the measure 77 above as \(r) + + tj_), 
where T7 + consists of the curves such that (cos 9, sin 9) ■ x 1 - > 0, and 77 _ consists of the curves such 
that (cos 9, sin 9) ■ ar 1 < 0, where x- 1 = (x2,— and the two flows 77 + and 77 _ can be shown to 
be still incompressible (see [4, Paragraph 4.1]). We will say more about these important examples 
later on, as more results on the theory will be available. 



3 The pressure field 

Bremer proved in [7] a surprising result: even though geodesies are not unique in general, given the 
initial and final conditions, there is a unique, up to an additive time-dependent constant, pressure 
field. The pressure field arises if one relaxes the incomprcssibility constraint, considering almost 
incompressible flows v. Denoting by p v the density produced by the flow, defined by 



(et)#i/ = p"(t, -)hd i.e. I 4>{u{t)) drj(uj) = §p v (t, •) dp D for all </> 



D 



we say that i> is almost incompressible if \\p v — lWc 1 < 1/2. 

Theorem 3.1 (Pressure as a Lagrange multiplier, [7, 1]) Let rj be optimal between 77 and 7. 
There exists a distribution p £ (C 1 )* such that 

st(v) + (p, p v - 1) > .5/(77) (5) 

for all almost incompressible flows v between 77 and 7 with p u (t,-) = 1 for t sufficiently close to 
and to 1. 

Using this result one can make first variations as follows: given a smooth field w(t, x), vanishing 
for t close to and 1, one can consider the family (X 1 ) of flow maps 

4-X*(e,x)=w(t,X t (e,x)), X t (0,x)=x 
de 

and perturb (smoothly) the paths w by u)(t) M> X (e, cj(t)) ~ oj(t) + e w(t,w(t)). Denoting by 

$ e : n*(D) -> n*(D), $ e (w,o)(t) := (X 4 (e, w(t)), a) , 

the induced perturbations in f2*(D), these in turn induce perturbations 77 e := ($ e )^77 of 77 which 
are almost incompressible. Then, the first variation gives 

f f (j(t) ■ ^j-w(t, uj(t)) dtdr}(u),a) + {p, divw) = 0. 
Jn*(D)Jo dt 

This equation uniquely determines Vp as a distribution, independently of the chosen minimizer 77: 
indeed, 77 enters in (5) only through stf(r\\ which obviously is independent of the chosen minimizer, 
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and so the above equation holds true for every minimizer r\. Since w is arbitrary, the first variation 
also leads to a weak formulation of Eulcr's equations 

d t vt(x) +div (v<8>v t (x)) + W x p(t,x) = 0, 

where v t and v <£> v t are implicitly defined by 

v t ^D = (e t )#(uj(t)ri), vWvt^D = (e*)#(w(i) ® u(t)rj). 

Observe that in general u <8> «t 7^ <8>«t- Indeed, since these models allow the passage of many fluid 
paths at the same point at the same time (i.e. branching and multiple velocities are possible), the 
product Vt (x) ® v t (x) of the mean velocity v t (x) with itself might be quite different from the mean 
value v v t (x) of the product. This gap precisely marks the difference between genuine distribu- 
tional solutions to Euler's equation and "generalized" ones (see also [4, Section 2 and Paragraph 4.4] 
for more comments on this fact). 



4 Necessary and sufficient optimality conditions 

In this section we study necessary and sufficient optimality conditions for Brcnier's variational prob- 
lem and its extensions. 

The basic remark is that any Borel integrable function q : [0, 1] X D — > R with j D q(t, •) d\iu = 
for every t £ [0,1] induces a null-lagrangian for the minimization problem, with the incompressibility 
constraint: indeed 

/ / q(t,u)(t)) dtdr/(uj,a) = / / q(t,x) dfiu(x) dt = 

Jn*(D) Jo Jo JD 

for any generalized incompressible flow r\. If we denote by 



the value function for the Lagrangian C q (ui) := J i|cO(£)| 2 — q(t,ui(t)) dt, we also have 

i i r 
-\uj(t)\ 2 -q(tMt))dtdri(u,a)> / ^/(a, h(a)) d» D (a) 

S2'(D) JO 1 JD 

for any incompressible flow rj between % and h. Moreover, equality holds if and only if r/-almost 
every (cj,o) is a c^ 1 -minimizing path. 

The following result, proved in [8, Section 3.6], shows that this lower bound is sharp with q = p, 
if p is sufficiently smooth. 

Theorem 4.1 Let u be a C 1 solution to the Euler equations in [0,T] x D, whose pressure field p 
satisfies 

(*) T 2 sup sup (Vlp(t,x)t;,0 <7T 2 . 

te[o,T] i£H, |£|<i 

Then the measure t] induced by u via the flow map is optimal on [0, T]. 

This follows by the fact that the integral paths of u satisfy iu(t) = —Wp(t,Lj(t)), and (*) implies 
that stationary paths for the action are also minimal for £ p . (This is a consequence of the one 
dimensional Poincare inequality \ii(t)\ 2 dt > ^ J Q T \u(t)\ 2 dt for all u : [0,T] — > R such that 
Jq udt = 0, see [5, Section 5] or [8, Proposition 3.2] for more details.) 
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The question investigated in [1] is: how far are these conditions from being necessary? C 1 
regularity or even one sided bounds on V 2 p are not realistic, so one has to look for necessary (and 
sufficient) conditions under much weaker regularity assumptions on p. 

From now on, we restrict for simplicity to the case D = T d . The following regularity result 
for the pressure field has been obtained in [2], improving the regularity Vp € M\ oc ((0, 1) x T d ) 
obtained in [8]. 

Theorem 4.2 For any 7, r\ E r(T <i ) the unique pressure field given by Theorem 3.1 belongs to 
Lf oc ((0,l);W(T d )). ' 

The above result says in particular that p is a function, and not just a distribution. This allows 
to define the value of p pointwise, which as we will see below will play a key role. 

In order to guess the right optimality conditions, we recall that the two main degrees of freedom 
in optimal transport problems are: 

• In moving mass from x to y, the path, or the family of paths, that should be followed; 

• The amount of mass that should be moved, on each such path, from x to y. 

The second degree of freedom is even more important in situations when more than one optimal 
path between x and y is available. As we will see, both things will depend on C p . But, since p is 
defined only up to negligible sets, the value of the Lagrangian C p on a path uj is not invariant in the 
Lebesgue equivalence class; furthermore, no local pointwise bounds on p are available (remember 
that p(t, •) is only a BV function, with BV norm in ij^O, 1)). Therefore, as done in [1], one has 
to: 

• Define a precise representative p in the Lebesgue equivalence class of p; it turns out that the 
"correct" definition is 

p(t,x) := liminf p(t, •) * </> e (a;), 

where p(t, •) * 4> e are suitable mollifications of p(t, •). Of course this definition depends on the choice 
of the mollifiers, but we prove that a suitable choice of them provides a well-behaved (in the sense 
stated in Theorem 4.3 below) function p. 

• Consider, in the minimization problem, only paths uj satisfying 

Mp(t,w(t)) e L\ oc {Q,l), 

where Mp(t, •) is a suitable maximal function of p(t,-) (see [1] for a more precise definition of the 
maximal operator). 

With these constraints one can talk of locally minimizing path uj for the Lagrangian Cp and, 
correspondingly, define a family of value functions 

cf : T d x T d -> [-00, +00], Mc(0,l), 

representing the cost of the minimal connection between x and y in the time interval (s,i): 

c/(x,y) := inf \\^{t)\ 2 -p(r,uj)dT : 

u(s) = x, u(t) = y, M p (t,uj{t)) e ^(M) j. 

With this notation, the following result proved in [1] provides necessary and sufficient optimality 
conditions. 

Theorem 4.3 Let r\ = r\ a ® /mt be an optimal incompressible flow between n = r\ a <g> pj and 7 = 
7a 65 Pt- Then 
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(i) r\ is concentrated on locally minimizing paths for Cp; 

(ii) for all intervals [s,t] C (0,T) ; for [if-a.e. a, the plan (e s , et)#r) a is c^'* -optimal, i.e. 

/ c/{x,y)d{e s ,e t ) # ri a < / c p \x,y)dX 

JT d xT d JT d xT d 

for any X G ^(T x T d ) having the same marginals of (e s , e t )#77 a . 

Conversely, if (i), (ii) hold with p replaced by some function q satisfying Mq G L\ oc ((0, 1); L 1 (Y d )) , 
then T] is optimal, and q is the pressure field. 

Notice that an optimal transport problem is trivial if either the initial or the final measure is a 
Dirac mass; therefore the second condition becomes meaningful when either (e s )#r) a or (et)#r] a is 
not a Dirac mass: it corresponds to the case when (e s , 7r a )^t?7 is not induced by a map, a phenomenon 
that cannot be ruled out, as we discussed in Section 2.2. In the example presented in Section 2.2 
the pressure field p{x) — \x\ 2 /2 is smooth and time-independent, but the initial and final conditions 
are chosen in such a way that a continuum of action-minimizing paths (4) between x and — x exists. 
As shown in [4], there are infinitely many incompressible flows connecting the identity map i to — i, 
which moreover induce infinitely many distributional solutions to the Euler equations [4, Paragraph 
4.4]. 

The results in [1] show a connection with the theory of action-minimizing measures, though in 
this case the Lagrangian J Q ^\uj(t)\ 2 — p(t,ui(t)) dt is possibly non-smooth and not given a priori, 
but generated by the variational problem itself. 

Here we see a nice variation on a classical theme of Calculus of Variations: a field of (smooth, 
nonintcrsccting) extremals gives rise both to minimizers and to an incompressible flow in phase 
space. Here, instead, we have a field of (possibly nonsmooth, or intersecting) minimizers which has 
to produce an incompressible flow in the state space. This structure seems to be rigid, and it might 
lead to new regularity results for the pressure field. 

Let us also recall that, as recently shown in [1 I], under a W^'P-regularity of the pressure p one 
can show that rj-a.e. uj solves the Euler- Lagrange equations and belongs to VF 2 ' p ([0, 1]) C C 1 ([0, 1]). 
This result is a first step towards the BV case, where one can still expect that the minimality of r) 
may allow to prove higher regularity on the minimizing curves u (like Co G BV). 
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